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We apply the perturbative chiral quark model (PCQM) to analyse low-energy nu- 
cleon properties: electromagnetic form factors, meson- nucleon sigma- terms and 
pion-nucleon scattering. Baryons are described as bound states of valence quarks 
surrounded by a cloud of Goldstone bosons (tt, K, n) as required by chiral symme- 
try. The model is based on the following guide lines: chiral symmetry constraints, 
fulfilment of low-energy theorems and proper treatment of sea-quarks, that is me- 
son cloud contributions. Analytic expressions for nucleon observables are obtained 
in terms of fundamental parameters of low-energy pion-nucleon physics (weak pion 
decay constant, axial nucleon coupling constant, strong pion-nucleon form factor) 
and of only one model parameter (radius of the nucleonic three-quark core). Our 
results are in good agreement with experimental data and results of other theoret- 
ical approaches. 



1 Introduction 

Hadron models set up to understand the structure of the nucleon should re- 
spect the constraints imposed by chiral symmetry. Spontaneous and explicit 
chiral symmetry breaking require the existence of the pion whose mass van- 
ishes in the limit of zero current quark mass. In turn, nucleon observables 
must receive contributions from pion loops. 

We recently suggested a a baryon model, the perturbative chiral quark 
model (PCQM), which includes relativistic quark wave functions and confine- 
ment and respects chiral symmetry. The PCQM was successfully applied to 
cr-term physijps U and extended teethe study of electromagnetic properties of 
the nucleon □ and ttN scattering u. Although the Lagrangian for this model 
fulfills the SU(2) x SU (2) current algebra, it is in general nontrivial to identify 
the explicit dynamics which leads to well-known predictions for various ob- 
servables. For instance, S-w&ve 7riV scattering at threshold is such a process 
which should be correctly described by a nucleon model involving the pion 
cloud. 
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2 Perturbative Chiral Quark Model 



2. 1 Effective Lagrangian and zeroth order properties 

Following considerations are based on the perturJpative chiral quari model 
(PCQM), a relativistic quark model suggested in □ and extended in El for the 
study of low-enetgy properties of baryons. Similar models have also been 
studied in Refs. EL The PCQM is based on the effective, chirally invariant 
Lagrangian C inv El 

C inv =4^- l°V(r) - S(r) 



U + E/t , U - E/t 
o + 7 5 S 



^Tr[^E/^E/t] (1) 



where ip is the quark field, U is the. .chiral field and F = 88 MeV is the pion 
decay constant in the chiral limit am. The quarks move in a self-consistent 
field, represented by scalar S(r) and vector V(r) components of an effective 
static potential providing confinement. The interaction of quarks, with Gold- 
stone bosons is introduced on the basis of the nonlinear cr-model We define 
the chiral field as U — exp[i(fr/F] using the exponentional parametrization of 
the nonlinear cr-model, where $ is the matrix of pseudoscalar mesons. When 
treating mesons fields as small fluctuations we do perturbation theory in the 
expansion parameter 1/F. Usually, the Lagrangian (|l]) is linearized with re- 
spect to the field $. Such an approximation is valid if we consider processes 
without free pions (e.g. nucleon mass shift due to the pion cloud) or with a 
single external pion field (e.g. pion-nucleon form factor) nr The resulting ap- 
proximate chiral invariance of the linearized Lagrangian Em guarantees both 
a conserved axial current (or PCAC in the presence of the meson mass term) 
and the Goldberger-Treiman relation. However, when we study ttN scattering 
the quadratic term in the expansion of the chiral field U should be kept. With 
the additional nonlinear term we reproduce □ the model-independent result 
for the 7riVjS'-wave scattering lengths as derived by derived by Weinberg and 
Tomozawa □. 

Here we concentrate on the application of the PCQM to the meson- 
nucleon cr-terms El and the electromagnetic properties of the nucleon H. We use 
the effective Lagrangian C e ff = + C-xSB including the linearized chiral- 
invariant term C^v an d a mass term C x sb which explicitly breaks chiral 
symmetry 

C l Z = # - S(r) - 7 °U(r)]^ + \(d^f - ^(r)z 7 5 |^, 
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C xSB = ~4>M^ - |Tr{$, {$, M}}, (2) 

where .M = diag{m, m, m s } is the mass matrix of current quarks (here and in 
the following we restrict to the isospin symmetry limit with m u = rrid = rh); 
B is the low-energy constant which measures the. vacuum expectation value 
of the scalar quark densities in the chkal limit E9. We rely on the standard 
picture of chiral symmetry breaking Ii3 and for the masses of pseudoscalar 
mesons we use the leading term in their chiral expansion (i.e. linear in the 
current quark mass): M% — 2rhB, M| = (m + m s )B, M% = |(m + 2m s )B. 
Meson masses satisfy the Gell-Mann-Oakes-Renner and the Gell-Mann-Okubo 
relation 3M 2 + M% — AM^. In the evaluation we use the following set of 
QCD parameters: m= 7 MeV, m s /m=25 and B = Af 2 + /(2m) = 1.4 GeV. 
Introduction of the electromagnetic field into the PCQM is accomplished 
by standard minimal substitution into Eq. (^). 

To describe the properties of baryons which are modelled as bound states 
of valence quarks surrounded by a meson cloud we formulate perturbation 
theory. In our approach the mass (energy) m™ re of the three-quark core 
of the nucleon is related to the single quark energy £o by m c ^ re = 3 • £o- 
For the unperturbed three-quark state we introduce the notation \<pQ > with 
the appropriate normalization < 4>q\4>q >= 1. The single quark ground state 
energy £ and wave function (w.f.) Uq(x) are obtained from the Dirac equation 

[-iaV + l3S(r) + V(r)-£o}uo(x) = 0. (3) 

The quark w.f. uq(x) belongs to the basis of potential eigenstates (including 
excited quark and antiquark solutions) used for expansion of the quark field 
operator ip(x). Here we restrict the expansion to the ground state contribution 
with ip(x) = boUo(x) exp(— iE^t), where bo is the corresponding single quark 
annihilation operator. In Eq. (||) the current quark mass is not included to 
simplify our calculational technique. Instead we consider the quark mass term 
as a small perturbation. 

For a given form of the potentials S(r) and V(r) the Dirac equation (|3|) 
can be solved numerically. Here, for the sake of simplicity, we use a variational 
Gaussian ansatz u for the quark wave function given by the analytical form: 



uq(x) = N exp 



2R 2 




(4) 



where N = [7r 3 / 2 i? 3 (l + 3p 2 /2)] 1 / 2 is a constant fixed by the normalization 
condition J d 3 x Uq(x) uq(x) = 1; \s, Xf, Xc are the spin, flavor and color 
quark wave functions, respectively. Our Gaussian ansatz contains two model 



NSTAR2001: submitted to World Scientific on February 1, 2008 3 



parameters: the dimensional parameter R and the dimensionless parameter p. 
The parameter p can be related to the axial coupling constant 17^4 calculated 
in zeroth-order (or 3q-core) approximation: 

5/ 2p 2 \ 5 1 + 2 7 

5A = 3l 1_ TTi7j = 3-^' (5) 



p- 



where 7 = 9(?a/10 — 1/2. The parameter i? can be physically understood as 
the mean radius of the three-quark core and is related to the charge, radius 
< r E > lo °f the P rc, ton in the leading-order (LO) approximation as B 



2 1 + |/0 2 V 2 

In our calculations we use the value <?yi=1.25 as obtained in ChPT i. We 
therefore have only one free parameter, R. In the numerical evaluation R is 
varied in the region from 0.55 fm to 0.65 fm. 

2.2 Perturbation theory and nucleon mass 

Following the Gcll-Mann and Low theorem we define the mass shift of the nu- 
cleonic three-quark ground state Amjv due to the interaction with Goldstone 
mesons as 

Am N = [ i6(t 1 )d 4 x 1 ...d 4 x n T[£ I (x 1 )...£ I (x n )]\</> >f (7) 

71=1 U ' J 

where Li = —S(r)ij) ij 5 ($/F) ip is the interaction Lagrangian between mesons 
and quarks treated as a perturbation and subscript " c" refers to contributions 
from connected graphs only. We evaluate Eq. (0) at one loop with o(l/F 2 ) 
using Wick's theorem and the appropriate propagators. For the quark field we 
use a Feynman propagator for a fermion in a binding potential. By restricting 
the summation over intermediate quark states to the ground state we get 

iG^(x,y) = < 4> \T{<tjj{x)iP(y)}\<Po > (8) 
-> u (x)u (y) exp[-i£ Q (x - y o )]0(x o - Vo)- 

For meson fields we use the free Feynman propagator for a boson field. Su- 
perscript "iV" in Eq. (^) indicates that the matrix elements are projected on 
the respective nucleon states. 

The total nucleon mass including one-loop corrections is given by 

m N =^£ Q + 1 m)+Y,d%Tl{Ml) ] d* N = ^,4 = ^ N ,dP N = ±<T N , (9) 
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where d% are the recoupling coefficients defining the partial contribution of 
the 7r, K and //-meson cloud to the mass shift of the nucleon. The contribution 
of a finite current quark mass ,1,0 the nucleon mass shift is taken into account 

le self-energj 

2 7 A 4 



5 rV 

perturbatively (for details see EJ). The self-energy operator H(M^) is given by 



where w^{p 2 ) — y/ Mg + p 2 is the meson energy with p = \p\ and F n j^^(p 2 ) 
is the ttNN form factor normalized to unity at zero recoil (p 2 = 0): 

W)-«p(-^){i + ^(i-^)}. (U, 

I?. 5 Meson nucleon sigma-terms and Feynman-Hellmann theorem 

The scalar density operators S i (i = u, d, s), relevant for the calculation 

of the meson-baryon sigma-terms in the PCQM, are defined as the partial 
derivatives of the model xSB Hamiltonian TC x sb = —jC x sb with respect to 
the current quark mass of i-th flavor mf. 



^PCQM _j_ 9TL x sb 



= Sr l + S? e a . (12) 



4 drrii 

Sy al is the set of valence-quark operators coinciding with the ones obtained 
from the QCD Hamiltonian. The set of sea-quark operators S^ ea arises from 
the pseudoscalar meson mass term a. To calculate meson-baryon sigma-terms 
we perform the perturbative expansion for the matrix element of the scalar 
density operator S PC( ^ M between unperturbed 3q-core states and then project 
it onto the respective baryon wave functions. 

As an example, the expression for a^N is given by 

*« N = m < p\S^ M + S P d CQM \ P> =3 7 m + J2 d%T(M 2 ) (13) 

where the first term of the right-hand side of Eq. ( |l3| ) corresponds to the 
valence quark, the second to the sea quark contribution. The vertex function 
r(M|) is related to the partial derivative of the self-energy operator II(M|) 
with respect to m: 

r(M|) = m^n(Mg). (14) 

Using Eqs. (|), (|l3| ) and ( |l4| ) we directly prove the Feynman-Hellmann the- 
orem o~ n N — ™ • dmN /dm. 
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3 Results 



3.1 Meson-nucleon sigma-terms 

We start our numerical analysis with the 7riV sigma-term. First, we restrict to 
the SU (2) flavor picture. Numerically, the contribution of the valence quarks 
to the irN sigma-term is 13.1 MeV, the contribution of sea quarks at order 
o(M%) is 66.9 ± 5.7 MeV. Here and in the following the error bars are due to 
variation of the range parameter R of the quark wave function (^) from 0.55 
to 0.65 fm. Taking into account higher-order contributions of the sea quarks 
we have the following result for the irN sigma-term a^ N (superscript it refers 
to the SU(2) flavor picture) of 

<jv = 43.3 ±4.4 MeV. (15) 

The contributions of kaon and 77-meson loops, a^ N and a^ N (superscripts 
K and r\ refer to the respective meson cloud contribution) are significantly 
suppressed relative to the pion cloud and to the valence quark contributions 
due to the energy denominators in the structure integrals. The same conclu- 
sion regarding the_suppression of K and ?7-meson loops was obtained in the 
cloudy bag model Ell. Numerically, kaon and 77-meson cloud contributions are 
<r% N = 1.7 ± 0.4 MeV and a^ N = 0.023 ± 0.006 MeV. For the irN sigma-term 
we have the following final value: 



^a% = 45 ±5 MeV. (16) 



Our result for the nN sigma-term is in perfect agreemcriLwith the value of 
a^N — 45 MeV deduced by Gasser, Leutwyler and Sainio E3 using dispersion- 
relation techniques and exploiting the chiral symmetry constraints. 

Next we discuss our prediction for the strangeness content of the nucleon 
Un which is defined in the PCQM as 

2< P \S™*"\p> 

JN < P \sr, cQM + s p d cQM \ P > [ ' 

The direct calculation of the strange-quark scalar density < p\Sf ® M \p > is 
completely consistent with the indirect one applying the Feynman-Hellmann 
theorem < p\Sf CQM \p >= dm N /dm s . The small value of y N = 0.076±0.012 
in our model is due to the suppressed contributions of kaon and 77-mesonJoops. 
Our prediction for is smaller than the value yN — 0.2 obtained in E3 from 
an analysis of experimental data on 7riV phase shifts. On the other hand, 
our prediction is quite close to the result obtained in the cloudy bag model 
y N w 0.05 Eil. 
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For the KN sigma-terms we obtain 

a u KN = 340 ± 37 MeV and a d KN = 284 ± 37 MeV, (18) 

which within uncertainties is consistent with values deduced in HBChPT and 
lattice QCD (see discussion in Ref. a). Hopefully, future DA$NE experiments 
at Frascati will allow for a determination of the KN sigma-terms and hence 
for a better knowledge of the strangeness content of the nucleon. 



3.2 Electromagnetic properties of the nucleon 

We extended our model analysis to the description of basic electromagnetic 
properties of the nucleon. Formal details can be found in Ref. □. We start 
with results for the magnetic moments of nucleons, /i p and fi n . For our set of 
parameters we obtain: 

Hp = 2.62 ±0.02, n n = -2.02 ±0.02, and — = -0.76 ± 0.01. (19) 

The leading order (LO, three-quark core) /J,p°, and next-to-leading order 
(NLO, meson cloud and finite current quark mass) contributions fJ.p LO , Hn LO 
to the magnetic moments are given by 

= 1.8 ±0.15, tf° = -lp%°, (20) 
^ LO = f i p - fj£° = 0.82 ± 0.13, ^ LO = Hn - = -0.82 ± 0.08. 
For the electromagnetic nucleon radii we obtain 

r p E = 0.84 ±0.05 fm, < r 2 >£= -0.036 ± 0.003 fm 2 , (21) 
r p M = 0.82 ± 0.02 fm, r n M = 0.85 ± 0.01 fm. 

The LO contributions to the charge radius of the proton (see Eq. (g)) and to 
the magnetic radii of proton and neutron are dominant 



r 



,p ;LO = 0.77 ±0.06 fm, r v ^° = r n ^ LO = 0.73 ± 0.06 fm. (22) 

For the neutron charge radius squared we get the observed (negative) sign, 
but its magnitude is smaller than the experimental value. As in the naive 
ST/ (6) quark model, the LO contribution to the neutron charge radius is zero 
and only one-loop diagrams give nontrivial contributions to this quantity: 

< r 2 > n E = -0.036 ± 0.003 fm 2 . (23) 

In Table 1 we summarize our results for the static electromagnetic proper- 
ties of the nucleon in comparison to experimental data. Results on (^e Q 2 
dependence of the electromagnetic form factors can be found in Ref. 
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Table l. Static nucleon properties. 



Quantity 


Our Approach 


Experiment 


M P 


2.62 ± 0.02 


2.793 




-2.02 ± 0.02 


-1.913 




-0.76 ± 0.01 


-0.68 


r E ( fm ) 


0.84 ± 0.05 


0.86 ± 0.01 


< r 2 >| (fm 2 ) 


-0.036 ± 0.003 


-0.119 ± 0.004 


r M ( fm ) 


0.82 ± 0.02 


0.86 ± 0.06 


rh (fm) 


0.85 ± 0.01 


0.88 ± 0.07 
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